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Abstract 



We point out that current experimental data for partial B — > irlv branching fractions 
reduce the theoretical input required for a precise extraction of \V u b\ to the form factor 
normalization at a single value of the pion energy. We show that the heavy-quark expan- 
sion provides a bound on the form factor shape that is orders of magnitude more stringent 
than conventional unitarity bounds. We find \V ub \ = (3.7±0.2±0.1) x [0.8/F+ (16 GeV 2 )]. 
The first error is from the experimental branching fractions, and the second is a conserva- 
tive bound on the residual form factor shape uncertainty, both of which will improve with 
additional data. Together with current and future lattice determinations of the form fac- 
tor normalization this result gives an accurate, model independent determination of \ V u b\. 
We further extract semileptonic shape observables such as |V^F + (0)| = 0.92±0.11±0.03 
and show how these observables can be used to test factorization and to determine low- 
energy parameters in hadronic B decays. 



1 Introduction 



Measuring the magnitude of the weak mixing matrix element V u b is important for constraining 
the unitarity triangle and testing the standard model of weak interactions. The exclusive 
determination of \V u b\ requires knowledge of the relevant heavy-to-light meson form factor and 
in the past this has led to significant model dependence in the result. First, the methods 
that were used to calculate the form factor, such as light-cone sum rules, quark models and 
quenched lattice calculations, all have unknown systematic errors. Second, each of these 
methods covers only part of the kinematic range; to obtain the total decay rate, the results 
were extrapolated using simplified parameterizations for the momentum dependence of the 
form factor. In the past year, the situation has improved dramatically: there are now several 
measurements of partial B — > ttIu branching fractions PjJ 121 El H] and the first results for 
the form factor from precision lattice simulations with dynamical light quarks have been 
presented jSJ E] . We show in this paper that if theoretical bounds on the form factor are taken 
into account, the experimental results for the partial branching fractions determine the shape 
of the form factor, to the point where no shape information at all is required from theory. This 
reduces the theoretical input for the determination of \V u b\ to a normalization of the relevant 
form factor, which can be taken at an energy within the range studied with current lattice 
simulations. For the first time, this allows an accurate, model independent, determination of 
\V u b\ from exclusive semileptonic decays. 

Bounds on the form factor can be derived via the computation of an appropriately chosen 
correlation function in perturbative QCD. By unitarity and analyticity, the resulting "disper- 
sive bound" constrains the behavior of the form factor in the semileptonic region |Hl El El] > 
and may be expressed as a condition on the coefficients in a convergent series expansion. While 
these bounded "series parameterizations" have been around for more than twenty-five years, 
many papers on the subject (in particular all experimental papers) have instead used simple 
pole forms to parameterize the form factor. In order to unify these descriptions, and to explain 
the dispersive bounds in a simple setting, we compare the class of series parameterizations 
emerging from the conventional dispersive bound analysis to the class of "pole parameteriza- 
tions" introduced in [TT) I12j. Both representations are exact, in the sense that the true form 
factor is guaranteed to be described arbitrarily well by a member of the class, and dispersive 
bounds can be established for both classes by power counting in the heavy-quark mass. With 
the resulting constraints in place, stable fits without truncation to a fixed number of parame- 
ters can be performed. In fact, for the series parameterization, we show that the bound given 
by simple heavy-quark power counting is orders of magnitude more stringent than the bound 
based on unitarity, thus providing much better control over the extraction of physical observ- 
ables from the data. We also derive three new exact sum rules for the coefficients appearing 
in the series representation of the form factor. 

The paper is organized as follows. In Section |21 we review the pole and series parameteriza- 
tions of the form factor. We list the experimental and lattice data to be used throughout the 
paper, and for later comparison we determine \V u b\ using three-parameter truncations of these 
parameterizations. Section El then introduces the bounds associated with each parameteri- 
zation, and establishes conservative estimates based on heavy-quark power counting for the 
bounded quantities. Section 0] examines the maximal precision for \V u b\ that can be reached 
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Figure 1: Experimental data for the partial B° — * tt + £~v branching ratios and fit result shown 
as solid line. The fit results from (J2J) with N = 1 and (jnj) with /c max = 2 are indistinguishable. 
Note that the experimental data is binned: PQ 12 E] give the result in three bins, while |4j gives 
the result in five g 2 -bins. We plot the value and error divided by the bin width at the average 
g 2 -value in each bin. For the three-bin results, we have slightly shifted the points to the left 
and right to increase visibility. 

with present data. We show that with a form factor determination at intermediate q 2 values, 
within the range studied in current lattice simulations, the experimental uncertainty on \V u i,\ 
is well below 10%, whereas a form factor determination near maximal q 2 would not translate 
into a precise value of \V u b\. We introduce three shape observables, |V^bF + (0)|, F+(0)/F + (0) 
and a, and discuss their sensitivity to the exact value of the bound. In Section EJ having 
established our procedure, we present final values for \V u b\, and for F + (0), in terms of a single 
lattice data point, _F + (16GeV 2 ). We extract the shape observables, which are determined by 
the experimental semileptonic data alone, and show how these observables provide important 
constraints on the factorization approach to hadronic B decays. 

2 Form factor parameterizations and extraction of \V u h\ 



Having restricted the shape of the q 2 spectrum, or equivalently, of the form factor, by exper- 
imental measurements, the central value and errors for \V u b\ are determined by varying the 
allowed form factor over all "reasonable" curves that are consistent with the data, and with a 
normalization of the form factor taken from theory at a given value (or multiple values) of q 2 . 
Defining this procedure precisely requires specifying a class of curves that contains the true 
form factor (to a precision compatible with the data), and that is sufficiently rich to describe 
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all variations impacting the observables under study. A statistical analysis along standard 
lines then determines central values and errors for the desired observable quantities. 

A starting point to isolate such a class of curves is the dispersive representation of the 
relevant form factor: 

F + (Q)/(1 - q) 1 r°° dt ImF + (t) 
1 - X- nj u t-q 2 -ie 



F + tf) = ' + ™* ; +r / dt^^-. (1) 



Here a is defined by the relative size of the contribution to i*+(0) from the 5* pole, and 
t ± = (m^i rn^) 2 . For massless leptons, the semileptonic region is given by < q 2 < i_. 
Equation (JTJ) states that, after removing the contribution of the B* pole lying below threshold, 
F + (q 2 ) is analytic outside of a cut in the complex g 2 -plane extending along the real axis from 
t + to oo, corresponding to the production region for states with the appropriate quantum 
numbers. 

One class of parameterizations keeps the B* pole explicit and approximates the remaining 
dispersion integral in ((TJ by a number of effective poles: 

F^) = F+[ * )l( \r a) + ±-^. ( 2) 

-L — — J — i 1 t — 3 — 

m-* B , k=l 7fe m^ t 

The true form factor can be approximated to any desired accuracy by introducing arbitrarily 
many, finely-spaced, effective poles. In the next section, we derive a bound on the magni- 
tudes, \pk\, of the coefficients of the effective poles. This allows a meaningful N — > oo limit, 
thus enabling us to investigate the behavior of the fits when arbitrarily many parameters are 
included. We find in actuality that current data cannot yet resolve more than one distinct 
effective pole in addition to the B* pole. Parameterizations of the above type are widely 
used to fit form factors. In particular, a simplified version of the N — 1 case, the so-called 
Becirevic-Kaidalov (BK) parameterization jTTj is used in many recent lattice calculations and 
experimental studies. As shown in ^2], this two-parameter form is overly restrictive since it 
enforces scaling relations which at small q 2 are broken by hard gluon exchange. The size of 
these hard-scattering terms, which appear at leading order in the heavy-quark expansion, is 
subject to some controversy and constraining their size is an important task. The parameter- 
ization of the form factors should allow for their presence. 

Another class of parameterizations is obtained by expanding the form factor in a series 
around some q 2 = to in the semileptonic region up to a fixed order, with the coefficients of this 
expansion as the fit parameters. The convergence of this simple expansion is very poor due 
to the presence of the nearby singularities at q 2 = m 2 B , and q 2 — t+. However, an improved 
series expansion of the form factor that converges in the entire cut g 2 -plane is obtained after 
a change of variables that maps this region onto the unit disc \z\ < 1. In terms of the new 
variable, F + has an expansion 

with real coefficients a*.. The variable z(q 2 ,to) maps the interval — oo < q 2 < t+ onto the line 
segment —1<z<1, with the free parameter to e (— oo, t + ) corresponding to the value of q 2 
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mapping onto z = 0. Points immediately above (below) the g 2 -cut are mapped onto the lower 
(upper) half-circle \z\ = 1. The function P(q 2 ) = z(q 2 ,m 2 B *) accounts for the pole in F + (q 2 ) 
at q 2 = m 2 Bt , while (p(q 2 ) is any function analytic outside of the cut. It is interesting to note 
that this reorganization succeeds in turning a large recoil parameter, (y ■ v ') max — 1 ~ 18, into 
a small expansion parameter. For example, for to — the variable z is negative throughout 
the semileptonic region and 



\ z \ = V> ■ ± 1 ' A w n 5 (4) 

max n < /— J V / 

^(v ■ v') max + 1 + V2 

where v and v' are the velocities of the B and 7r mesons. The same size, but for positive z 
is obtained for t = £_. By choosing the intermediate value t = ^+(1 — v/1 — the 
expansion parameter can be made as small as |z| max ~ 0.3. A second class of parameterizations 
is obtained by a truncation of the above series: 

F+{q2) = P(^(^E^W (5) 



0; 



fc=0 



As discussed in the next section, it is conventional to take 
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l^j ^fr#j (T^#j (t^- C) 



With this choice, a bound a 2 . < 1 is obtained by perturbative methods. 1 Together with 
the restriction \z\ < 1, this allows a meaningful k max — > oo limit. In actuality, we find that 
current data can only resolve the first three terms in the series (|5|i. 

Figure n shows the available experimental data on the partial branching fraction dT(B° — > 
ix + i~ v) I dq 2 . The CLEO £Q, Belle [2] and BaBar [3] collaborations have measured this branch- 
ing fraction in three separate g 2 -bins and BaBar [3] has presented a measurement using five 
g 2 -bins. The correlation matrix is included in our fits for the data in pQ. For the remaining 
data, g 2 -bins are taken as uncorrelated. In order to extract \V u b\ we also need the normaliza- 
tion of the form factor. The B — > 7r vector form factors have recently been determined by the 
Fermilab Lattice and by the HPQCD [Oj collaborations in lattice simulations with dynam- 
ical fermions. The preliminary results of these calculations give F + (16 GeV 2 ) = 0.81 ±0.11 5J 
and F + (16GeV 2 ) = 0.73 ± 0.10 0. 2 Although the lattice calculations give the form factor 

x For different to ? the expansion parameters, z = z(t, to) and z' = z(t, t' ), and expansion coefficients, 
ak = ak(to) and a' k = <Zfc(td), are related by the Mobius transformation: 

, z(t ,t' ) + z 



1 + z(t ,t' n )z 



ft is easily verified that the sum of squares of coefficients is invariant under such a transformation, af, = 
J2k °fc 2 ' as guaranteed by the construction of tf>, see @. 

2 The parameterization J2J) with N = 1 has been used to interpolate to the common g 2 -point, and for 
definiteness the errors are taken as those from the nearest points: q 2 = 1 5. 87 GeV 2 0, with statistical and 
systematic errors added in quadrature, and q 2 = 1 6. 28 GeV 2 [B]. 
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Figure 2: 68% (dark) and 95% (light) confidence limits for \V u b\ determined by fitting the 
parameterizations (0) or (0) to experimental data in pQ, [2|, |S] and jl], with the single lattice 
data point -F + (16 GeV 2 ) = 0.8 ± 0.1. Results from (j2J) and (0) are indistinguishable. The plot 
on the right shows \V u b\ for fixed .F + (16GeV 2 ) = 0.8 as a function of the relative uncertainty 
on the form factor. 

at several different g 2 -values, the correlations between different points are not available and 
it is difficult to quantify the uncertainty on the shape. Anticipating the analysis of Section 01 
where we determine the range of q 2 that best exploits the experimental shape information, 
we use F + (16GeV 2 ) = 0.8 ±0.1 as our default value for the form factor normalization. Note 
that we have avoided any theoretical biases concerning the form factor shape. Performing a 
X 2 fit yields \V u b\ = 3.7^05 x 10~ 3 for both the parameterization (J2J) with N — 1, and (0) with 
k — 9 

""max 

Figure El shows the 68% and 95% confidence limits for \V u b\ as a function of the value 
and uncertainty of the form factor at q 2 = 16 GeV 2 . The form factor normalization is the 
dominant error in the determination of \V u b\; if the quantity F + (16GeV 2 ) would be known 
exactly, the uncertainty on \V u b\ would drop to approximately 6%. The quality of the fit is 
equally good for both parameterizations, with x 2 — 12.0. 3 The extracted value of \V u b\ is 
insensitive to the choice of the free parameter t - Setting <p{q 2 ) — 1 in (JSJ) also has negligible 
impact, and similarly adding more lattice input points does not substantially change the result 
if the dominant lattice errors are correlated. The effect of allowing additional terms in the 
parameterizations (0) and © is investigated in the following sections. We will find that the 
result for the value and uncertainty of \V u b\ from the simple parameterizations used in this 

3 Note that all three-bin measurements determine the same observable quantities. The minimal x 2 obtained 
from the three-bin measurements is 5.0 for 9 — 3 degrees of freedom. This value measures the (good) agreement 
between the three-bin measurements, and should be subtracted from the total in order to obtain a measure of 
agreement between the data and the parameterizations. The resulting quality of our fit is good: 12.0 — 5.0 for 
9 — 4 degrees of freedom. 
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section is not appreciably altered if additional terms are included. 



3 Form factor bounds 

To make a fully rigorous determination of \V u i,\, the truncation to the three-parameter classes 
of curves considered in the previous section requires justification. For instance, if the neglected 
terms in (J2J) or (jSJ) conspired to produce a sharp peak in the form factor at precisely the value 
of the lattice input point, then the integrated rate would be overestimated, and the value of 
\V u b\ underestimated. To prevent this from happening requires some bound on the perversity 
of allowed form factor shapes. In practice, we would like to ensure that our extraction of 
physical observables is "model- independent" by allowing for arbitrarily many parameters, i.e., 
taking N — » oo in (J2J and A; max — > oo in (jSJ). Retaining predictive power then demands that a 
bound be enforced on the parameters appearing in these expansions. 

To bound the coefficients pk in the expansion (j2J, we introduce a decomposition of the 
integration region, t+ < t± < • ■ ■ < tjv+i < oo, and define 

1 f th+1 dt T „ . , t k 
Pk = - — ImF+(i) , 7^ = —2- . (7) 

Since F + (t) ~ t^ 1 at large t, it follows that 

5>*l<-/ -\F + (t)\=R, (8) 

and this is the desired bound. The integral in (JSJ) is dominated by states with t — t + ~ m b A, 
where F + ~ m^" 2 ', so that the quantity R is parametrically of order (A/m?,) 1 / 2 , with A a 
hadronic scale. To be sure that the bound deserves the model-independent moniker, one 
should use a very conservative estimate. In our fits we will use R < vTo and R < 10, i.e., we 
allow for an addition factor of 100 or 1000 beyond the dimensional estimate R 2 ~ A/rrib ~ 0.1. 

The coefficients in the expansion (jSJ) can be bounded by requiring that the production 
rate of Bti states, described by the analytically continued form factor, does not overwhelm 
the production rate of all states coupling to the current of interest (in this case, the vector 
current uj^b). The latter rate is computable in perturbative QCD using the operator product 
expansion (for a pedagogical discussion, see e.g. [TU] ). The function in (jHJ) was chosen such 
that the fractional contribution of Bn states to this rate is given at leading order by 

f>l = ^ f -\Hz)P(z)F + (z)\ 2 = ^ r^[(t - t + )(t - t^\F + (t)\ 2 = A. (9) 
fc=o 1X1 J z Jt + 

In the heavy quark limit, the leading contributions to the integral A in @ are of order (A/m?,) 3 
and arise from two regions: the region close to threshold, t — t + ~ m^A, where the pion has 
energy E ~ A and the form factor scales as F + ~ m^ 2 ; and the region t — t + ~ m 2 , where 

E ~ rrib and F + ~ m b 3 ^ 2 (for a discussion of the form factor scalings, see ^2])- The region of 
very high energies t>mj, where F + ~ 1/t, gives a sub leading contribution. 
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Since, by definition, the fraction is smaller than unity, it is conventional to take the loose 
bound A < 1, which does not make use of scaling behavior in the heavy-quark limit. Clearly 
this bound leaves much room for improvement; from its scaling behavior, we expect A to be on 
the order of a few permille. This implies that higher-order perturbative and power corrections 
in the operator product analysis introduce negligible error, as noticed in [T3]. It is also easy 
to see that the dispersive bounds by themselves do not impose tight constraints on the form 
factor shape. Since the scale of the coefficients is set by ao ~ m^ 3 ^ 2 , even with the optimal 
choice |z| m ax ~ 0.3, the dispersive bounds allow the relative size of higher-order terms in the 
series, \akZ k /ao\, to be of order unity up to k « 4, and to contribute significantly for even 
higher k. This situation for heavy-to-light decays such as B — > 7r contrasts with that for 
heavy- heavy decays such as B — ► D [TH IT5] , where the bound is parametrically of order unity 
(counting m c ~ m;, ^> A). For this case, the scale of the coefficients is set by ao ~ m", and 
with |^| m ax ~ 0.06 the bound ensures that only the first few terms in the series are required 
for percent accuracy. 

To put the dispersive bounds in perspective, it may be useful to emphasize that establishing 
an order-of- magnitude bound on any integral of the form f dt k(t)\F + (t)\ 2 for some k(t) would 
yield an equally valid, bounded, parameterization, with a new <p(t) constructed from k(t) as 
in @ and 0. Similarly, bounded pole parameterizations (j2J) are obtained by establishing an 
order-of-magnitude bound on any integral of the form J^dt k(t)\F + (t)\ for some kit), as in 
(J7J) and ©. Focusing attention on the special case of © and is justified only to the extent 
that the bound © is sufficiently restrictive, and to the extent that similar or tighter bounds 
cannot be conservatively estimated by other means. 

It is interesting to note that the two bounds are not equivalent. The bound Ylk \Pk\ ^ 
R < oo uses the fact that the asymptotic form factor can be evaluated in perturbation theory, 
where the scaling F + (t) ~ t" 1 is found at large t. This condition is not automatically satisfied 
by the series parameterization (jSJ), which as seen from © requires only F+(t) < t 1/2 at large 
t. Imposing the proper large-t behavior yields the sum rules 

= 0, n = 0,1,2. (10) 

Z->1 

To our knowledge, the above sum rules have not been discussed in the literature. On the other 
hand, all pole parameterizations "violate" the bound Y^T=o a 2 . = A < oo for the simple reason 
that the integral in @ is not well defined for these parameterizations, because F + (t) has poles 
on the integration contour. 

The bounds discussed here are associated with the behavior of the form factor above thresh- 
old. Since we are interested in the form factor in the semileptonic region, these higher-energy 
properties are useful only to the extent that they can help to constrain the form factor in this 
region. Incorrect high-energy behavior therefore does not imply that a given parameterization 
cannot be used to describe low-energy data. For instance, the effective poles in (J2J) could be 
smeared into finite- width effective resonances in order to make the integral in converge; 
however, the semileptonic data is very insensitive to such fine-grained detail, and this modifi- 
cation has a very minor impact on the fits. Similarly, unless the bound © is close to being 
saturated, the coefficients a k for moderately large k in the series parameterization © can 
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Figure 3: A% 2 = 1 region for \ V u j,\ for an infinitely precise form factor determination at a single 
g 2 -value. The plot assumes that the form factor yields the central value \V u b\ = 3.7 x 10~ 3 . 
The darkest band is obtained for J2 k al < 0.01, while the two lighter bands correspond to 

J2 k 4 < o- 1 and Ek a l < L 

be tuned to satisfy the sum rules (|TU|) . or equivalently, to make the integral in (jHJ) converge. 
However, the semileptonic data becomes insensitive to terms z k for large k, and again such a 
modification has little impact on the fits. Thus, while at some level the bound (jHJ) will con- 
strain the parameters in the series parameterization (JSJ), and the bound Q will constrain the 
parameters in the pole parameterization (j2J), we restrict attention to the constraints imposed 
by (jSJ) on the pole parameterization, and by (jUJ) on the series parameterization. 

4 Parameterization uncertainty and shape observables 

With the bounds in place, it is straightforward to generalize the fits in Section El to include 
arbitrarily many parameters. Imposing the very conservative bound \p k \ < 10, we observe 
that additional poles in the class of parameterizations (j2J) have essentially no impact on the 
central value and errors for \V u b\. Similarly, using the very conservative bound a\ < 1 in 
((HJ), we find that the inclusion of higher order terms beyond fc max = 2 has negligible impact 
on \V u b\. The errors are dominated by the lattice input point, and both the central value and 
errors are not changed significantly from the iV = 1 or /c max = 2 fits in Section El 

In order to isolate the uncertainty on the form factor shape inherent to the present data, 
we show in Figure El the minimum attainable error on \V u b\, assuming exact knowledge of the 
form factor at one g 2 -value. Results are shown for the parameterization (JSJ), using various 
bounds J2k a k < 0.01, 0.1 and 1. As the figure illustrates, points in the intermediate range of 
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Table 1: Fit results for form factor shape parameters using the pole parameterization (j2J). 

g 2 lead to the smallest uncertainty on \V u b\, and for these points, the \V u b\ extraction is not 
very sensitive to even the order of magnitude of the chosen bound, with the minimum error 
varying from approximately 6% to approximately 8% as the bound is relaxed from 0.01 to 1. 
It should be noted that a better understanding of correlations in the experimental data would 
be necessary when probing this level of precision. The curves in Figure E2 are also indicative 
of the impact of additional theory inputs. Performing the fits with data points at different 
g 2 -values in addition to the default F + (16GeV 2 ) shows that a point at q 2 = would require 
< 10% error to significantly decrease the error on \V u b\, while even exact knowledge of the 
form factor at q 2 = £_ has almost no impact. 

In the remainder of this section we consider observables which are more sensitive to the 
shape of the form factor and investigate the role played by the bounds in these cases. In par- 
ticular, we extract the form factor and its first derivative at q 2 = 0, as well as the residue at the 
B* pole, which is directly related to the parameter a, as in These quantities are interest- 
ing in their own right. The form factor at zero momentum transfer, normalized as |\4fe|-f+(0), 
is an important input for the evaluation of factorization theorems for charmless two-body 
decays such as B — > tttt. The derivative of the form factor at q 2 = 0, conveniently normalized 
as {m 2 B — m 2 )F^_(0)/F + (0), determines the quantity 5 measuring the ratio of hard-scattering 
to soft-overlap terms in the form factor Finally, the value of (1 — a)" 1 is proportional 
to the coupling constant gB*Bn- The observable quantities |14&|F + (0), F|(0)/F + (0) and a 
are independent of the form factor normalization, and hence are determined solely by the 
experimental data. 

In TablesnandlU we show how the results for the shape observables change when additional 
parameters are included. For the pole parameterization (J2J) we perform fits with N = 1,2 and 
3 poles in addition to the B* pole. (The case N = 1 was studied in (T^j.) To help stabilize 
the fits, we impose a minimum spacing of the poles 'jk+i ~ Ik > and a maximum pole 

position, 7& < N + 1. For the polynomial parameterization JHJ), we set A; max = 2, 3 and 4. We 
perform each of the fits with two different bounds — a very loose model-independent bound, 
and a more stringent bound that relies on the scaling behavior of the bounded quantity in the 
heavy-quark limit. Given a value of the bound, a central value and errors are determined by 
taking the limit of large iV in (J2J), or large fc max in (jSJ). The sequence converges once the size 
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Table 2: Fit results for form factor shape parameters using the series parameterization (JHJ) 
with t = 0. 



of the neglected terms is constrained by the bound to lie below the sensitivity of the chosen 
observable. 

The quantities |Ki&|-F+(0), F|(0)/F + (0) and a exhibit different sensitivities to the bounds. 
This is to be expected, since sharp bends in the fitted curve at the endpoints allowed by the 
additional terms can have strong effects on the slope, or on the residue of the B* pole, but 
are not constrained tightly by the data. Imposing only very loose bounds therefore leads to 
large uncertainties for these quantities. 

It is instructive to examine the relation between observables and expansion coefficients. 
At to = the quantities /(0), a, (3 and 5 studied in [TJ] are related to the coefficients by 

1/2 



/(0) = F + (0) 



16ao 



[l+m^f 2 1 + 



1 + /T 



B 



m; 



dF+ 



F+(0) dq 2 



a 




—ai 1 — m„ 
4a 1 + m-x ' 4 1 + 

(l + m, + A) 2 (l + y/^:) 3 ' x 
4(l+m T ) 2 (A + 2 



tub + m n ) 2 — m 2 B *, and hats denote quantities in units of m#. The heavy- 
quark scaling laws for f(0) and 1 + — 5 are special cases of the general law ~ m^ 3 ^ 2 , 

~ fn^^ 2 when 

scaling violations are neglected. Similarly, the scaling law for (1 — a) -1 
behavior A -1 / 2 ~ ml J 4 for the sum appearing in the last line of (jll|) . 



where A 2 
quark sea 

obtained by taking k derivatives in (j3J), and noticing that d n F + /d(q 2 ) n 



q 2 =0 

translates into the 



5 Results and discussion 

In order to extract the most precise value of \V U b\, it is important to make full use of the 
existing experimental data for B — > nlu that determines the form factor shape. To emphasize 
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this point, the analysis was done here using no shape information at all from theory, but only 
a normalization at one g 2 -point. Our results make it clear that the limiting factor in the 
determination of \V u b\ is currently the form factor normalization, with very small uncertainty 
associated with the form factor shape. Similar conclusions are implicit in other recent works. 
For example, in the reduction in error compared to methods employing only total exper- 
imental branching fractions is due almost entirely to the inclusion of shape information from 
experiment, and not to the inclusion of additional theory input points. In experimental 
data is combined with simple parameterizations of the form factor shape to constrain the 
hadronic input parameters appearing in sum rule estimates of the form factor. In contrast to 
these and other previous works, we have avoided any theoretical biases concerning the form 
factor shape. 

In practical terms, the parameterizations (j2J, with N — 1, and with fc ma x = 2, are 
sufficient for describing the current generation of semileptonic data, in the sense that the 
addition of more parameters does not significantly improve the fits. To provide rigorous 
error estimates it is necessary to allow for arbitrarily many additional parameters within the 
dispersive bounds (jHJ) and 0. For "global" quantities like \V u b\ it is possible to show by 
imposing only the very loose bounds ^2 k \pk\ < 10 in (J2J), or J2k a l < 1 in © that the 
extracted values are actually insensitive to the addition of more parameters. With a single 
lattice input value F + (16 GeV 2 ) = 0.8 ± 0.1, we find 

iKcfcl = 3.7 ± 0.2 ±g;« ± 0.1 

F+(0) = 0.25 ± 0.04 ± 0.03 ± 0.01 

The first error is experimental, the second is theoretical from the lattice input, and the third 
is due to the uncertainty in the form factor shape. For definiteness, the central values in 
flT^j) are obtained using the parameterization (JSJ with ^2 k a\ < 0.01, and the third error is 
very conservatively estimated by adding the maximum variation of the boundaries of the la 
interval induced by relaxing the bound to J2k a t < 1- 

For less global quantities, like the slope of the form factor at q 2 = 0, the very loose bounds 
(JHJ) and (jnj) are not sufficient to tightly constrain the impact of arbitrarily many additional 
parameters. In this case we adopt more realistic estimates for the bounds, and find 

10 3 |Kfe^+(0)| = 0.92 ± 0.11 ± 0.03 , 

1.5 ±0.6 ±0.4, (13) 
8 1* ± 7 . 

The first error is experimental, and the second is due to uncertainty in the form factor shape 
(these quantities are independent of the form factor normalization). The central values in 
()13|) are again obtained using the parameterization (jSJ) with a\ < 0.01, and the shape 
error is conservatively estimated by adding the maximum variation of the boundaries of the 
la interval when the bound is relaxed to Ylik a \ < 0-1- 



(3.7 ±0.2) x 



0.8 

F+(16GeV : 



2\ ' 



(0.25 ±0.04) x 



F + (16GeV 2 
08 



(12) 



, n(Q) 

*+(0) 



a 
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While the conventional dispersive bound approach provides an elegant means of demon- 
strating formal convergence properties with the minimal assumption of form factor analyticity 
and the convergence of an operator product expansion, some caution is required in order to 
avoid misinterpreting the results. Firstly, for certain observables, e.g. \V u b\, the fits are much 
more tightly constrained by the data than by the unitarity-based dispersive bound. This leads 
to the happy conclusion that the errors on \ V u b\ do not depend on the chosen parameterization 
or the exact value of the bound, and the analysis lends itself to a straightforward statistical 
interpretation. Secondly, other important observables, such as the slope of the form factor, 
are sensitive to the addition of more parameters than can be constrained by the data, but are 
allowed by the unitarity bound. Since this bound is overestimated, presumably by orders of 
magnitude, a reliance on this procedure would lead to the pessimistic conclusion that almost 
no information at all can be extracted from the data for these quantities. In such cases, we 
propose to use tighter bounds, which follow from the scaling behavior of the bounded quantity 
in the heavy quark limit. 

Apart from establishing order-of-magnitude estimates for the bounds in (jHj) and by 
heavy-quark power counting, none of the above analysis relies on heavy-quark, large-recoil 
or chiral expansions, or on the associated heavy-quark, soft-collinear or chiral effective field 
theories. However, the semileptonic data can be used to test predictions from these effective 
field theories, and to determine low-energy parameters that can be used as inputs to the 
calculation of other processes. For example, using the experimental result Br(5~ — > 7r~7r°) = 
(5.5 ± 0.6) x 10~ 6 Hl| together with |K 6 |F+(0) from (JT3J), we find 



T(B- 



dT(B° -> n+l~v)/dq 2 \ 



,2- 



:0 



0.76 ±°; 22 ± 0.05 GeV 2 , (14) 



where the first error is experimental, and the second is due to the form factor shape uncertainty 
in (|13|). Such ratios provide a strong test of factorization JS]- The leading-order prediction for 
this ratio, corresponding to the "naive" factorization picture where hard-scattering corrections 
are neglected, yields 167r7 2 |Kd| 2 (Ci + C 2 ) 2 /3 = 0.62 ± 0.07 GeV 2 . This uncertainty includes 
only the effects of varying the renormalization scale of the leading-order weak-interaction 
coefficients (HI] between m^jl and 2m;>. This may be compared to the prediction of Beneke 
and Neubert j2H| who use QCD factorization theorems for two-body decays to work beyond 
leading order and include the effects of hard-scattering terms, obtaining for the same ratio, 
0.66 

iaos GeV 2 . The uncertainty in their prediction is dominated by the uncertainty in the 
light-cone distribution amplitudes (LCDAs) of the B- and 7r-mesons. Bauer et al. j^U EI| 
evaluate the same factorization theorems using a different strategy: they use experimental 
results for other B — > tttt decays to determine the part involving the LCDAs from data, which 
is possible if all power corrections, and perturbative corrections of order a s (m&), are neglected. 
For the ratio (JHJ) they find 1.27^^ 

GeV 2 , where we display only experimental errors. The 
semileptonic data provides important information on otherwise poorly constrained hadronic 
parameters entering these processes. 

As a second application, the parameter S measuring the relative size of hard-scattering and 
soft-overlap contributions in the B — > tt form factor can be related to the slope of the form 
factor at q 2 = [12 . Extrapolated to zero recoil, the lattice calculations in [SJ E] give for the 
slope of the F form factor, (3 = [(m| - m„) 2 F^O) / F+iO)]- 1 = 1.2 ± 0.1. Together with (JT3J) 
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this yields 



m 2 B — m 2 ( dF+ 



F + (0) V dq 2 



dq 2 



0.4 ±0.6 ±0.1 ±0.4, (15) 



where the first error is experimental, the second is theoretical from the lattice determination of 
(3, and the third is due to the form factor shape uncertainty in Qlflj) . Establishing the relative 
size of the hard-scattering and soft-overlap contributions from the semileptonic data provides 
another important input to factorization analyses of hadronic B decays. The above result for 
S does not unambiguously establish 5^0 which signals the presence of hard-scattering terms, 
but it disfavors the opposite scenario, 5 ~ 2, where the form factor is completely dominated 
by hard-scattering. More data will help reduce both the experimental and shape-uncertainty 
errors for this quantity. 

As a third application, the form factor F + (0) and shape observable a determine the cou- 
pling constant Qb*b-k via 

—————— = = 2.0 ± 0.2 ± 1.7 , (16) 

2m b* 1 — a 

where the first error is experimental, the second is theoretical from the lattice form factor 
normalization, and the final error is due to the form factor shape uncertainty, determined as 
in (|13|). Since the semileptonic data is concentrated at small q 2 , it is not very sensitive to the 
detailed structure of the sub-threshold pole and dispersive integral in ((TJ). In fact, the data 
do not yet definitively resolve a distinct contribution of the B* pole, although the opposite 
scenario — dominance by the B* pole in (JTJ) — is ruled out [T2| . 

Our implementation of the bounds in (JSJ) and Q could be formalized in terms of stan- 
dard methods of constrained curve fitting [22] ■ In this language, we have enforced a "prior" 
probability function which is constant if the parameters obey the bound on \pk\ or a 2 ,, 
and zero otherwise. For simplicity, we then performed a x 2 fit, assuming sufficient statistics 
that the data is Gaussian distributed. The resulting error estimates should be conservative. 
Firstly, this prior allows equal probability for parameter values that are near the bound, even 
though we believe such values are increasingly unlikely. Other prior functions may be con- 
sidered — for example, in the case of the series parameterization (jSJ), a Gaussian prior on 
the variable (— log 10 J2 k a 2 .), with mean and standard deviation of order unity. Secondly, in 
estimating errors based on A% 2 , we neglect the fact that bounds enforce restrictions that 
renormalize the probability distributions, and to the extent that the bounds are relevant, this 
tends to overestimate errors. As a simple example, if an absolute bound happened to coin- 
cide with the boundary of the "la" interval obtained for an observable based on A% 2 = 1, 
we would estimate that the observable was within the interval with only ~ 68% confidence, 
whereas the bounds guarantee this with 100% confidence. In a more refined analysis, a direct 
evaluation of the statistical integrals could account for such boundary effects. An alternative 
procedure employed in and generalized in [23J to include shape information from experi- 
ment, has a slightly more complicated statistical interpretation. Here theory information on 
the form factor, combined with the dispersive bounds, is used to generate a statistical sample 
of "envelopes", each consisting of the curves defined at each g 2 -point by the minimum and 
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maximum values that the form factor can take. (Note that some curves may be ruled out by 
the bounds, yet allowed by the envelopes, which are generated by extremizing point-by-point 
in q 2 .) This sample of envelopes is then combined with experimental branching fractions to 
determine a distribution for \V u b\ or other observables. Working in terms of parameters 
allows the experimental and lattice data to be treated on the same footing, and yields a more 
straightforward interpretation of the constraints enforced by the bounds. Fortunately, these 
complications play an extremely minor role in the case of \V u b\. As illustrated by Figure ^ the 
errors are very nearly Gaussian, and nearly identical results are obtained using different pa- 
rameterizations, and widely different values for the bounds. A more refined statistical analysis 
might be useful for those shape observables that show sensitivity to the bounds, to extract as 
much information as possible from the experimental data. 

In summary, we have shown that the form factor shape information necessary for a pre- 
cise extraction of \V u b\ is now entirely determined from experiment. Rather than relying on 
theoretical models for this shape, the current and future experimental data can instead be 
used as a precision tool for testing theory predictions and determining hadronic parameters 
in other processes. For example, the ratio in (jl4j) should be predicted with good accuracy 
from the factorization approach to hadronic B decays, and can be even more firmly estab- 
lished once the hard-scattering contribution in (|15|) is determined more precisely from data. 
The methodology employed here in B decays can be validated in the analogous situation of 
semileptonic D decays, where experiment and lattice cover the entire range of q 2 . Note that 
we only used lattice input for the form factor at a single g 2 -value, to avoid theoretical biases 
on the form factor shape, and to emphasize the conclusion that the shape is determined by 
experiment; however, studying the form factor shape provides an important test of lattice 
calculations. Our results show that with improved lattice data, an exclusive measurement of 
| V u b | that rivals or even surpasses the inclusive determination is possible. 
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